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We show that the transmission through single and double (5-function potential barriers of strength 
P = VWb/hvp in bilayer graphene is periodic in P with period n. For a certain range of P values 
we find states that are bound to the potential barrier and that run along the potential barrier. 
Similar periodic behaviour is found for the conductance. The spectrum of a periodic succession of 
5-function barriers (Kronig-Penney model) in bilayer graphene is periodic in P with period 2tt. For 
P smaller than a critical value Pc, the spectrum exhibits two Dirac points while for P larger than Pc 
an energy gap opens. These results are extended to the case of a superlattice of (5-function barriers 
with P alternating in sign between successive barriers; the corresponding spectrum is periodic in P 
with period tt. 
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I. INTRODUCTION 

Graphene, a one-atom thick layer of carbon atoms, has 
become a research attraction pole since its experimental 
discovery in 2004^. Since carriers in graphene behave like 
relativistic and chiral massless fermions with a linear-in- 
wave vector spectrum, many interesting features could be 
tested with this material such as the Klein paradox^''^, 
which was recently observec3, the anomalous quantum 
Hall effect, etc., see Ref. "Ffor two recent reviews. The 
effort to realise this Klein tunnelling through a potential 
barrier also lead to other interesting features, such as 
resonant tunnelling through double barrier^. With the 
possibility to fabricate devices with single-layer graphene, 
bilayer graphene has also been extensively investigated 
and been shown to possess extraordinary electronic be- 
haviour, such as a gapless spectrum, in the absence of 
bias, and chiral carrier^^. Many of these nanostructures 
could be given another functionality if based on bilayer 
instead of single-layer graphene. 

The electronic band structure can be modified by the 
application of a periodic potential and/or magnetic bar- 
riers. Such superlattices (SLs) are commonly used to 
alter the band structure of nanomaterials. In single- 
layer graphene already a number of papers relate their 
work to the theoretical understanding of such periodic 
structure^^"^. Much less experimental and theoretical 
work has been done on bilayer graphene^* 

We will study the spectrum, the transmission, and the 
conductance of bilayer graphene through an array of po- 
tential barriers using a simple model: the Kronig-Penney 
(KP) model^^, i.e. a one-dimensional periodic succession 
of (S-function barriers on bilayer graphene. The advantage 
of such a model system is that, 1) a lot can be done ana- 
lytically, 2) the system is clearly defined, 3) and it is pos- 
sible to show a number of exact relations. The present re- 
search is also motivated by our recent findings for single- 
layer graphengi^, where very interesting and unexpected 
properties were found, for instance, that the transmis- 



sion and energy spectrum are periodic in the strength of 
the (5-function barriers. Surprisingly, we find that for bi- 
layer graphene similar, but different, properties are found 
as function of the strength of the (5-function potential 
barriers. Due to the different electronic spectra close to 
the Dirac point, i.e., linear for graphene and quadratic 
for bilayer graphene, we find very different transmission 
probabilities through a finite number of barriers and very 
different energy spectrum, for a superlattice of (5-function 
barriers, between single-layer and bilayer graphene. 

The paper is organised as follows. In Sec. II we briefly 
present the formalism. In Sec. Ill we give results for the 
transmission and conductance through a single (5-function 
barrier. We dedicate Sec. IV to bound states of a single 
(5-function barrier and Sec. V to those of two such barri- 
ers. In Sec. VI we present the spectrum for the KP model 
and in Sec. VII that for an extended KP model by con- 
sidering two (5-function barriers with opposite strength in 
the unit cell. Finally, in Sec. VIII we make a summary 
and concluding remarks. 



II. BASIC FORMALISM 

We describe the electronic structure of an infinitely 
large flat graphene bilayer by the continuum nearest- 
neighbour, tight-binding model and consider solutions 
with energy and wave vector close to the K {K') point. 
The corresponding four-band Hamiltonian and eigen- 
states ^ are 
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with t: = px + ipy {px,y — ~ifidx,y) and p the mo- 
mentum operator. We apply one-dimensional potentials 
V{x, y) = V{x) and consequently the wave function can 
be written as '4>{x, y) = '4){x)e'^^^'^ with the momentum in 
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the y-direction a constant of motion. Solving the time- 
independent Schrodinger equation H^p — we obtain, 
for constant V{x, y) = V, the spectrum and the eigen- 
states. The latter are given by Eq. (A5) in App. [A| and 
the spectrum by Eq. (|A2| 



e = u + 1/2 ± v/l/4 + fc2, 
e = u- l/2± Vl/4 + fc2, 



(2) 



where we used the dimensionless variables, e = E/t±, 
u — V/t_\_, X xt±/hvF, ky — ^ hvFky/t± and e' — s — u; 
vf = 10^ m/s, and t± = 0.39 eV expresses the coupling 
between the two layers. 

For later purposes we also give the frequently used two- 
band Hamiltonian 



n 



Vf 



TT^' 
t2 



and the corresponding spectrum 

E-v = ±{vF^hyt^){ki 



(3) 
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As seen, there are qualitative differences between the two 
spectra (compare Eq. Q with Eq. that will be re- 
flected in those for the transmission and conductance in 
some of the cases studied. As the approximation of the 
two-band Hamiltonian is only valid for E <^ t±^, we can 
expect a qualitative difference with the four-band Hamil- 
tonian if \E\ > i_L. 



III. TRANSMISSION THROUGH A 
(5-FUNCTION BARRIER 

We assume jE' — y| < t± outside the barrier such that 
we obtain one pair of localised and one pair of travelling 
eigenstates in the well regions characterised by wave vec- 
tors a and /3, where a is real and /3 imaginary, see App. 
El Consider an incident wave with wave vector a from 
the left (normalised to unity); part of it will be reflected, 
with amplitude r, and part of it will be transmitted with 
amplitude t. Then the transmission is T = \t\'^. Also, 
there are growing and decaying evanescent states near 
the barrier, with coefficients eg and e^, respectively. The 
relation between the coefficients can be written in the 
form 
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This leads to a system of linear equations that can be 
written in matrix form 
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with Nij the coefficients of the transfer matrix 



noting the matrix in Eq. ([6| 
coefficients from (t, r, e^, 
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we can evaluate the 

T 



by Q 

r, ed, eg,)^ - Q-^ (l, 0, 0, O) 
As a result, to obtain the transmission ampli- 
tude t it is sufficient to find the matrix element 



V33J 



We model a i5-function barrier as the limiting case of 
a square barrier, with height V and width Wj, shown in 
Fig. [1] represented by V{x) = VQ{x)Q{Wi, - x). The 
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FIG. 1: Schematics of the potential V(x) of a single square 
barrier. 

transfer matrix N for this (5-function barrier is calculated 
in appendix [B] and the limits V 00 and Wb — >■ are 
taken such that P — VWb/hvF is kept constant. 

The transmission T — |ip for a real and /3 imaginary 
is obtained from the inverse amplitude. 
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, (7) 



where /i = {e+1/2) /a and ly = {e—l/2)//3. Contour plots 
of the transmission T are shown in Figs. |2ja) and |2jb) 
for strengths P = 0.257r and P = O.TStt, respectively. 

The transmission remains invariant under the trans- 
formations: 



1) P^ P + rnr, 



2) ky 



kii 



(8) 



The first property is in contrast with what is obtained in 
Ref. 3 In the latter work it was found, by using the 2x2 
Hamiltonian, that the transmission T should be zero for 
ky K, Q and E <Vq, while we can see here that for certain 
strengths P — mr there is perfect transmission. The last 
property is due to the fact that ky only appears squared 
in the expression for the transmission. Notice that in 
contrast to single-layer graphene the transmission for e « 
is practically zero. The cone for nonzero transmission 
shifts to £ = 1/2(1 — cosP) with increasing P till P — 
TT. An area with T = appears when a is imaginary, 
i.e., for e'^ + e — ky < (as no propagating states are 
available in this area, we expect bound states to appear) . 
From Figs. [2]^a), 2j^b) it is apparent that the transmission 
in the forward direction, i.e., for ky « 0, is in general 
smaller than 1; accordingly, there is no Klein tunneling. 
However, for P = mr, with n an integer, the barrier 
becomes perfectly transparent. 
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For P = riTT we have V = liVF{mT/Wb). If the elec- 
tron wave vector is fc = rm/Wb its energy equals the 
height of the potential barrier and consequently there is 
a quasi-bound state and thus a resonancd^. The condi- 
tion on the wave vector implies Wb = nX/2 where A is 
the wavelength. This is the standard condition for Fabry- 
Perot resonances. Notice though that the invariance of 
the transmission under the change P ^ P + mr is not 
equivalent to the Fabry-Perot resonance condition. 
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FIG. 2: (Color online) Contour plot of the transmission for 
P = 0.257r in (a) and P = 0.757r in (b). In (b) the bound 
state, shown by the red curve, is at positive energy. The white 
area shows the part where a is imaginary. The probability 
distribution |^/)(a;)p of the bound state is plotted in (c) for 
various values of ky and in (d) for different values of P. 



From the transmission we can calculate the conduc- 
tance G given, at zero temperature, by 



G/ Go 



Tr/2 



(9) 
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where Go = {4,e^ /2T:h)[El, + t^EFY'"^ /hvp] the angle 
of incidence (j) is determined by tani^ = ky/a. It is not 
possible to obtain the conductance analytically, therefore 
we evaluate this integral numerically. 

The conductance is a periodic function of P (since the 
transmission is) with period tt. Fig. |3] shows a contour 
plot of the conductance for one period. As seen, the con- 
ductance has a sharp minimum at £ = 1/2(1— cos P): this 
is due to the cone feature in the transmission which shifts 
to higher energies with increasing P. Such a sharp min- 
imum was not present in the conductance of single-layer 
graphene when applied to the same (5-function potential 
barriei^^'. 



FIG. 3: (Color online) (a) Contour plot of the conductance 
G. (b) Slices of G along constant P. 



IV. BOUND STATES OF A SINGLE 
(5-FUNCTION BARRIER 

The bound states here are states that are localised in 
the x-direction close to the barrier but are free to move 
along the barrier, i.e. in the y-direction. Such bound 
states are characterised by the fact that the wave function 
decreases exponentially in the x direction, i.e., the wave 
vectors a and /3 are imaginary. This leads to 





(10) 



which we can write as 



= Q 




(11) 



where the matrix Q is the same as in Eq. ^. In or- 
der for this homogeneous algebraic set of equations to 
have a nontrivial solution, the determinant of Q must be 
zero. This gives rise to a transcendental equation for the 
dispersion relation 



det Q = iVii7V33 - TViaiVai = 0, 
which can be written explicitly 

(a - pfe 



(12) 



[cos P+iiJ, sin P] [cos P+iv sin P] + - 



sin^ P = 0, 
(13) 



This expression is invariant under the transformations 
P + n-TT, 

> -ky, 



1) P- 

2) ky 

3) (e,P)^ (-£,^-P), 



(14) 



Furthermore, there is one bound state for fcy > and 
7r/2 < P < TT. For P < Tr/2 we can see that there is 
also a single bound state for negative energies from the 
third property above. From this transcendental formula 
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one can find the solution for the energy e as function of 
ky numerically. We show the bound state by the solid 



red curve in Fig. [2|(b). This state is bound to the po- 
tential in the x direction but moves as a free particle 
in the y direction. We have two such states, one that 
moves along the +y direction and one along the ~y di- 
rection. The numerical solution approximates the curve 
e = cosP [-1/2-1- (1/4 -t- fc2)i/2]. If one uses the 2 x 2 
Hamiltonian one obtains the dispersion relation given in 
Appendix [C] by Eq. ( |C1[ ). By solving this equation one 
finds for each value of P two bound states one for pos- 
itive and one for negative ky. Moreover, for positive P 
these bands have a hole like behaviour and for negative 
P an electron like behaviour. Only for small P do these 
results coincide with those from the 4x4 Hamiltonian. 

The wave function ip{x) of such a bound state is char- 
acterised by the coefficients e^i, 6^2 on the left, and edi 
and 6(12 on the right side of the barrier. We can obtain 
the latter coefficients by using Eq. (11), by assuming 
Bgi — 1 and afterwards normalising the total probability 
to unity in dimensionless units. The wave function ^{x) 
to the left and right of the barrier can be determined from 
these coefficients by using Appendix [A| In Figs. |2jc), (d) 
we show the probability distributioiiJi/;(a::)P of a bound 
state for a single ^-function barrier: in (c) we show it 
for several ky values and in (d) for different values of P. 
One can see that the bound state is localised around the 
barrier and is less smeared out with increasing ky. No- 
tice that the bound state is more strongly confined for 
P ~ tt/2 and that |?/;(x)p is invariant under the trans- 
formation P — > TT — P. 



V. TRANSMISSION THROUGH TWO 
<5-FUNCTION BARRIERS 

We consider a system of two barriers, separated by 
a distance L, with strengths Pi and P2, respectively, as 
shown schematically in Fig. |4] We have L — ^ Lt±/hvF = 
0.59261L/nm which for i = 10 nm, vp = 10^ m/s, and 
t± — 0.39 eV equals 5.9261 in dimensionless units. The 
wave functions in the different regions are related as fol- 
lows 



V(x) 

(1) (2) t (3) 



where S' = GM{1)G ^ represents a shift from x=0 to 
x=L and the matrices 5*1 and S2 are equal to the matrix 
N' of Eq. ([B3| with P = Pi and P = P2, respectively. 
Using the transfer matrix Af = Q^^SiS'S2GM{L) we 
obtain the transmission T = 

In Fig. 5] the transmission T{e, ky) is shown for par- 
allel (a), (bj and anti-parallel (b), (c) (S-function barriers 
with equal strength, i.e., for |Pi| — IP2I, that are sepa- 
rated by L = 10 nm, with P — 0.257r in (a) and P = O.Stt 
in (b). For P — n/2, the transmission amplitude t for 
parallel barriers equals —t for anti-parallel ones and the 
transmission T is the same, as well the formula for the 
bound states. Hence panel (b) is the same for parallel 
and anti-parallel barriers. The contour plot of the trans- 
mission has a very particular structure which is very dif- 
ferent from the single-barrier case. There are two bound 
states for each sign of ky, which are shown in panel (d) for 
parallel and panel (e) for anti-parallel barriers. For anti- 
parallel barriers these states have mirror-symmetry with 
respect to e = but for parallel barriers this symmetry is 
absent. For parallel barriers the change P — >■ tt — P will 
ffip the spectrum of the bound state. The spectrum of the 
bound states extends into the low-energy transmission 
region and gives rise to a pronounced resonance. Notice 
that for certain P values (Figs. |5ja) and[5](d)) the energy 
dispersion for the bound state has a camelback shape for 
small ky, indicating free propagating states along the y 
direction with velocity opposite to that for larger ky val- 



ues. Contrasting Fig. |2jb) with Fig. |5[d)-(e)we see that 
the free-particle like spectrum of Fig.T2[b) for the bound 
states of a single (5-function barrier is strongly modified 
when two (5- function barriers are present. 

From the transfer matrix we find that the transmission 
is invariant under the change P — ^ P + mr and ky — )> —ky 
for parallel barriers, which was also the case of a single 
barrier, cf. Eq. (14 1. In addition, it is also invariant, for 



anti-parallel barriers, under the change 

P ^ TT - P. 



(16) 



The conductance G is calculated numerically as in the 
case of a single barrier. We show it for (anti-)parallel 
(5-function barriers of equal strengths in Fig. [6] The 
symmetry G(P + rnr) — G(P) of the single barrier 
conductance holds here as well. Further, we see that 
for anti-parallel barriers G has the additional symmetry 
G(P) — G{tt — P) as the transmission does. 



VI. KRONIG-PENNEY MODEL 



FIG. 4: A system of two (5-function barriers with strengths 
Pi and P2 placed a distance L apart. 



We consider an infinite sequence of equidistant S- 
function potential barriers, i.e., a superlattice (SL), with 
potential 



V{x) = P^(5(x-nL). 



(17) 



MO) - SiMo), Mo) = s'ML), 

M^) ^ S2ML), Mo) = SiS'S2ML), (15) 



As this potential is periodic the wave function should be 
a Bloch function. Further, we know how to relate the 
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FIG. 5: (Color online) Panels (a), (b), and (c): contour plots of the transmission through two J-function barriers of equal 
strength P = |Pi| = IP2I separated by a distance L — 10 nm. For parallel barriers we took P = 0.257r in (a) and P — O.Stt in 
(b). For anti-parallel barriers results are given for P — 0.5n in (b) and P = 0.257r in (c). The solid red curves in the white 
background region is the spectrum for the bound states. Panels (d) and (e) show the dispersion relation of the bound states 
for various strengths \P\, respectively, for parallel and anti-parallel barriers. The thin black curves delimit the region where 
bound states are possible. 
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FIG. 6: (Color online) Contour plot of the conductance of 
two 5-function barriers with strength IP2I = |Pi| ~ P and 
inter-barrier distance L — 10 nm. Panel (a) is for parallel 
barriers and panel (c) for anti-parallel barriers. Panels (b) 
and (d) show the conductance, along constant P, extracted 
from panels (a) and (c), respectively. 



coefficients Ai of the wave function before the barrier 
with those (^3) after it, see Appendix [b] The result is 



V'(O), Ai^J^As, 



(18) 



with kx the Bloch wave vector. From these boundary 
conditions we can extract the relation 



M{L)A3^UA3 



(19) 



with the matrix A4{x) given by Eq. (A4|. The determi 



nant of the coefficients in Eq. ( |19| must be zero, i.e., 

det[e-'''''^M{L) - = 0. (20) 



If ky = 0, which corresponds to the pure ID case, one 
can easily obtain the dispersion relation because the first 
two and the last two components of the wave function 
decouple. Two transcendental equations are found 

cos kxL = cosQ!icosP+ - f — h —\ sinaLsinP, (21a) 

2 V £ a/ 

cos kxL = cos PL cos P + 2 ~^ ^) ^'^ ^'^^ ^' ^'^^^^ 
Since /3 is imaginary for < £^ < tj_, we can write Eq. 

^L^sinh|/3|LsinP, (22) 
2|P|£ 



(pl|3) as 

cos fcri = cosh \B\L cos P 



which makes it easier to compare with the spectrum of 
the KP model obtained from the 2x2 Hamiltonian, see 
Eq. (3). The latter is given by the two relations 

cos k^L — cos kL -f (P/2k) sin kL, (23a) 
cos k^L = cosh — (P/2k) sinh kL, (23b) 

with K — yfe. This dispersion relation, which has the 
same form as the one for standard electrons, is not 'peri- 
odic in P and the difference from that of the four-band 
Hamiltonian is due to the fact that the former is not valid 
for high potential barriers. One can also contrast the dis- 
persion relations (21) and (|23| with the corresponding 



one on single-layer graphene 



cos k^L — cos XL cos P + sin XL sin P, 



(24) 



where A = E/{hvF). This dispersion relation is also pe- 
riodic in P. 

In Fig. |7]we plot slices of the energy spectrum for ky = 
0. There is a qualitative difference, between the four- 
band and the two-band approximation for P = ir. Only 
when P is small does the difference between the two ID 
dispersion relations become small. Therefore, we will no 
longer present results from the 2x2 Hamiltonian though 
it has been used frequently due to its simplicity. The 
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present results indicate that one should be very careful 
when using the 2x2 Hamiltonian in bilayer graphene. 

Notice that for P = 0.257r the electron and hole bands 
overlap and cross each other close to \ky\ « 0.5(7r/L). 
That is, this is the spectrum of a semi-metal. These 
crossing points move to the edge of the Brillouin zone 
(BZ) for P = TT resulting in a zero-gap semiconductor. 
At the edge of the BZ the spectrum is parabolic for low 
energies. 

For ky 7^ 0, the dispersion relation can be written ex- 
plicitly in the form 
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FIG. 7: (Color online) Slices of the spectrum of a KP SL 
with L = 10 nm along k^, for ky — 0, with P = 0.257r in 
(a) and (b) and P = tt in (c) and (d). The results in (a) 
and (c) are obtained from the four-band Hamiltonian (1) and 
those in (b) and (d) from the two-band one (3). The solid and 
dashed curves originate, respectively, from Eqs. ( 21 1 23 1) and 
(l2Tb[23t). 



cos2k^L + CiCosk^L + Co/2 ^ 0, (25) 

where 

Ci = -2(cos aL + cos /3i) cos P - {da + dp) sin P, (26) 
and 

Co = {2 + kl / e'^) + {2 - k^/s'^) cos aL cos PL 



lis' 



-e^{2e^ - l)]sinaLsin/3L/2a/3e^ 



- ^kl/e^ - {2 + kl/e^) cos aL cos 13 L 

+ [2e^ -l/2-kl + k^/e^] sinaLsin/3L/a/3| cos2P 
-I- [da COS /3L + dp cos aL] sin 2P, (27) 

with da — {2e + 1) sin aL/a and dp — {2e — 1) sin (3L//3. 
The wave vectors a = [e^ + e — ky]^^^ and /3 = [e^ — e — 
fc^]^/^ are pure real or imaginary. If /3 becomes imaginary, 
the dispersion relation is still real (/3 — i\/3\ and sin /3L — 



jsinh |/3|J^). Further, if a becomes imaginary, that is for 
a — > i\a\, the dispersion relation is real. The dispersion 
relation has the following invariance properties: 

1) E{k,,ky,P) = E{k,,ky,P + 2mr), (28a) 

2) E{k,,ky,P) = -E{T:/L-k,,ky,TT-P),{2%h) 

3) E{k,,ky,P) = E{k,,~ky,P). (28c) 

In Fig. |8] we show the lowest conduction and highest va- 

(b) 




FIG. 8: (Color online) SL Spectrum for L — IQ nm, the 
lowest conduction and highest valence band for P — 0.257r in 
(a) and P — 0.5tv in (b), are shown. 
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FIG. 9: (Color online) SL spectrum for L = 10 nm. The 
dashed blue, solid red, and dash-dotted purple curves are, 
respectively, for strengths P — O.lvr, P — 0.257r, and P = 
O.Stt. (a) shows the spectrum vs k^ for ky = while (b) 
shows it vs ky for k^ at the value where the bands cross. The 
position of the touching points and the size of the energy gap 
are shown in (c) as a function of P. The dash-dotted, blue 
curve and the solid, black curve show kx,o and the energy 
value of the touching points, respectively. For P > P^ a. gap 
appears; the energies of the conduction band minimum and of 
the valence band maximum are shown by the red and purple, 
solid curve, respectively. 

lence bands of the energy spectrum of the KP model for 
P = 0.257r in (a) and P = O.Stt in (b). The former has 



7 



two touching points which can also be viewed as over- 
lapping conduction and valence bands as in a semi-metal 
and the latter exhibits an energy gap. In Fig. |9] slices of 
Figs. [Sj^a), (b) are plotted for ky — 0. The spectrum of 
bilayer graphene has a single touching point at the origin. 
When the strength P is small, this point shifts away from 
zero energy along the axis with ky = and splits into 
two points. It is interesting to know when and where 
these touching points emerge. To find out we observe 



that at the crossing point both relations (21 1 should be 
fulfilled. If these two relations are subtracted we obtain 
the transcendental equation 

= (cosaL-cos/3L)cosP+(l/2) (d„ - dp)smP, (29) 

where = (2e + 1) s'maL/a and = (2£— 1) sin/3i//3. 
We can solve Eq. (28) numerically for the energy e. For 
small P and small L this energy can be approximated by 
e = P/L. Afterwards we can put this solution back into 
one of the dispersion relations to obtain k^- 

In Figs, ^a), (b) we show slices along the kx axis 
for ky — and along the ky axis for the kx value of a 
touching point, kxfi- We see that as the touching points 
move away from the K point, the cross sections show a 
more linear behaviour in the ky direction. The position of 
the touching points is plotted in Fig. ^c) as a function 
of P. The dash-dotted blue curve corresponds to the 
value of kx,o (right y axis), while the energy value of the 
touching point is given by the black solid curve. This 
touching point moves to the edge of the BZ which occurs 
for P = Pc ~ 0.4257r. At this point a gap opens (the 
energies of the top of the valence band and of the bottom 
of the conduction band are shown by the lower purple and 
upper red solid curve, respectiv ely) and increases with P. 
Because of property 2) in Eq. 28 b) we plot the results 
only for P < 7r/2. We draw attention to the fact that 
the dispersion relation differs to large extent for large P 
from the one that results from the 2x2 Hamiltonian 
given in Appendix [C] This is already apparent from the 
fact that the dispersion relation does not exhibit any of 
the periodic in P behaviours given by Eqs. (28i) and 

An important question is whether the above period- 
icities in P still remain approximately valid outside the 
range of validity of the KP model. To assess that we 
briefly look at a square-barrier SL with barriers of fi- 
nite width Wb and compare the spectra with those of the 
KP model. We assume the height of the barrier to be 
V/hvF = P/Wb, such that VWb/hvF = P- The SL pe- 
riod we use is 50 nm and the width Wb = 0.05L = 2.5 nm. 
For P = tt/2 the corresponding height is then V ^ t±. 
To fit in the continuum model we require that the poten- 
tial barriers be smooth over the carbon-carbon distance 
which is a « 0.14 nm. In Fig. [TOjwe show the spectra for 
the KP model and this SL. Comparing (a) and (b) we 
see that for P small the difference between both models 
is rather small. If we take P — tt/2 though, this dif- 
ference becomes large, especially for the first conduction 
and valence minibands, as shown in panels (c) and (d). 



The latter energy bands are fiat for large ky in the KP 
model, while they diverge from the horizontal line (E=0) 
for a finite barrier width. From panel (f), which shows 
the discrepancy of the SL minibands between the exact 
ones and those obtained from the KP model, we see that 
the spectra with P = G.27r are closer to the KP model 
than those for P = 7r/2. Fig.[To{^e) demonstrates that the 
periodicity of the spectrum in P within the KP model, 
i.e., its invariance under the change P ^ P + 2m:, is 
present only as a rough approximation away from it. 



VII. EXTENDED KRONIG-PENNEY MODEL 

In this model we replace the single J-function barrier 
in the unit cell by two barriers with strengths Pi and 
—P2- Then the SL potential is given by 

V{x) = Pi'Y^5{x-nL)-P2'Y^5{x-{n+l/2)L). (30) 

n n 

Here we will restrict ourselves to the important case of 
Pi — P2. For this potential we can also use Eq. (20 1 
of Sec. |VI| with the transfer matrix J\f replaced by the 
appropriate one of Sec. [Vj 

First, let us consider the spectrum along ky = which 
is determined by the transcendental equations 



cos kxL — cos aL cos^ P + Da sin^ P, (31a) 
cos kxL — cos PL cos^ P + Dp sin^ P, (3 lb) 



with Dy ~ [(7^ + £^) COS7L - 



7^ -f-e' 



/47^ 



It is 



more convenient to look at the crossing points because 
the spectrum is symmetric around zero energy. This fol- 
lows from the form of the potential (its spatial average 
is zero) or from the dispersion relation (31a): the change 
e — )■ —e entails a ^ (5 and the crossings in the spectrum 
are easily obtained by taking the limit e — > in one of 
the dispersion relations. This gives the value of kx at the 
crossings 



kxfi = ± arccos[l - (L^/8) sin^ P]/L, 



(32) 



and the crossing points are at (e, kx, ky) — (0, ±kxfi, 0). 
If the kxfi value is not real, then there is no solution at 
zero energy and a gap arises in the spectrum. From Eq. 
(31a I we see that for sin^ P > 16/ a band gap arises. 

In Fig. [11] we show the lowest conduction and highest 
valence band for (a) P = 0.1257r, and (b) P = 0.257r. If 
we make the correspondence with the KP model of Sec. 
V we see that this model leads to qualitatively similar 
(but not identical) spectra shown in Figs, [sja) andjsj^b): 
one should take P twice as large in the corresponding 
KP model of Sec. |V]in order to have a similar spectrum. 
Here we have the interesting property that the spectrum 
exhibits mirror symmetry with respect to e = which 
makes the analysis of the touching points and of the gap 
easier. 

12 we plot the kx value (dash-dotted, blue 
le touching points kxfl versus P, if there is 
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FIG. 10: (Color online) Spectrum of a SL with L = 50 nm, (a) and (b) are for P — 0.2n and (c) and (d) are for P = n/2. (a), 
(c) and (e) are for a rectangular-barrier SL with Wb = 0.05L and u = P/Wb, while (b) and (d) are for the KP model, (e) shows 
the spectrum for u corresponding to P = (1/2 + 2)n; the dashed curves show the contours of the spectrum in (c) for P = n/2. 
(f) Shows the discrepancy of the SL minibands between the exact ones and those obtained from the KP model, averaged over 
k space (where we used kyL/n = 6 as a. cut-off). The conduction (valence) minibands are numbered with positive (negative) 
integers. 




FIG. 11: (Color online) The first conduction and valence 
minibands for the extended KP model for L = 10 nm with 
P = 0.1257r in (a) and P = 0.257r in (b). 



0.5 




P/JI 



FIG. 12: (Color online) Plot of the ik^^fi values, for which the 
minibands touch each other, as a function of P (dash-dotted, 
blue curve), and the size of the band gap Egap (solid, red 
curve). The calculation is done for the extended KP model 
with L = 10 nm. 



no gap, and the size of the gap Egap (solid, red curve) 
if there is one. The touching points move toward the 
BZ boundary with increasing P. Beyond the P value for 
which the boundary is reached, a gap appears between 
the conduction and valence minibands. 



VIII. CONCLUSIONS 

We investigated the transmission through single and 
double (5-function potential barriers on bilayer graphene 
using the four-band Hamiltonian. The transmission and 
conductance are found to be periodic functions of the 
strength of the barriers P = VWb/hvp with period tt. 
The same periodicity was previously obtained for such 
barriers on single-layer graphene^-. We emphasise that 
the periodicity obtained here implies that the transmis- 
sion satisfies the relation T{kx, ky, P) = T{kx^ ky, P+m:) 
for arbitrary values of kx, ky, P, and integer n. In 
previous theo retical work on graphen^^Sl and bilayer 
graphenJ^mil Fabry-Perot resonances were studied and 
T — 1 was found for particular values of a, the elec- 
tron momentum inside the barrier along the x axis. For 
a rectangular barrier of width W and Schrodinger-type 
electrons, Fabry-Perot resonances occur for aW = mr 
and E > Vq as well as in the case of a quantum well 
for E > 0, Vq < 0. In graphene, because of Klein tun- 
nelling, the latter condition on energy is not needed. Be- 
cause a depends on the energy and the potential barrier 
height in the combination E — Vq, any periodicity of T 
in the energy is equivalent to a periodicity in Vq if no 
approximations are made, e.g., E <C Vq, etc. Although 
this may appear similar to the periodicity in P, there 
are fundamental differences. As shown in Ref. 21, the 
Fabry-Perot resonances are not exactly described by the 
condition aW — mr (see Fig. 3 in Ref. 21) while the pe- 
riodicity of T in the effective barrier strength is exactly 
riTT. Furthermore, the Fabry-Perot resonances are found 
for T = 1 , while the periodicity of T in P is valid for any 
value of T between and 1 . 

Further, we studied the spectrum of the KP model 
and found it to be periodic in the strength P with period 
27r. In the extended KP model this period reduces to 
TT. This difference is a consequence of the fact that for 
the extended SL the unit cell contains two ^-function 
barriers. These periodicities are identical to the one 
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found earlier in the (extended) KP model on single-layer 
graphene. We found that the SL conduction and valence 
minibands touch each other at two points or that there 
is a energy gap between them. In addition, we found a 
simple relation describing the position of these touching 
points. None of these periodic behaviours results from 
the two-band Hamiltonian; this clearly indicates that the 
two-band Hamiltonian is an incorrect description of the 
KP model in bilayer graphene. In general, results de- 
rived from these two tight-binding Hamiltonians agree 
well only for small energie^. The precise energy ranges 
are not explicitly known and may depend on the partic- 
ular property studied. For the range pertaining to the 
four-band Hamiltonian ab-initio result^^H indicate that 
it is approximately from —1 eV to + 0.6 eV. 

The question arises whether the above periodicities in 
P survive when the potential barriers have a finite width. 
To assess that we briefly investigated the spectrum of a 
rectangular SL potential with thin barriers and compared 
it with that in the KP limit. We showed with some ex- 
amples that for specific SL parameters the KP model is 
acceptable in a narrow range of P and only as a rough ap- 
proximation away from this range. The same conclusion 
holds for the periodicity of the KP model. 

The main differences between the results of this work 
and those of our previous one, Ref. TT", follows. In 

contrast to monolayer graphene we found here that: 
1) The conductance for a single (5-function potential bar- 
rier depends on the Fermi energy and drops almost to 
zero for certain values of E and P. 2) The KP model 
(and its extended version) in bilayer graphene can open 
a band gap; if there is no such gap, two touching points 
appear in the spectrum instead of one. 3) The Dirac line 
found in the extended KP model in single-layer graphene 
is not found in bilayer graphene. 



be obtained, for constant V{x, y) = V, by progressive 



elimination of the unknowns in Eq. (All and solution 
of the resulting second-order differential equations. The 
result for the spectrum is 



e = u+l/2± Vl/4 + fc2, 
e = u-l/2± ^l/4 + fc2. 



(A2) 



The unnormalised eigenstates are given by the columns 
of the matrix QA4, where 



g = 



( 1 


1 


1 


M 


fa 
J + 


fa 






-1 


-1 


1 




V/- 


fa 
J + 




-n) 



(A3) 



with f^''^ = -i{ky±i{a,P))/e'; a = [e'^ +e' - fc2]i/2 and 



(3 = 
by 



e' — kyY^^ are the wave vectors. Ai is given 



M 



f^iax 




V 










\ 





—i/Sx J 



(A4) 



The wave function in a region of constant potential is a 
linear combination of the eigenstates and can be written 




(A5) 



We can reduce its complexity by the linear transforma- 
tion ^{x) -J> TZ'li{x) where 
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7^ = 




(A6) 



which transforms \E'(a;) to '^{x) = (1/2)(V'a — V's'i V's ~ 
ijjA' , ipA + '4'B', "0-8 + V'A' )'^- Then the basis functions are 
given by the columns of QA4 with 



Appendix A: Eigenvalues and eigenstates for a 
constant potential 

Starting with the Hamiltonian ([l]) for a one- 
dimensional potential V{x,y) — V(x), the time- 
independent Schrodinger equation Hijj = Eijj leads to 



-i{drc - ky)^B = e'i^A - ipB', 
-i{dx + ky)ipA = s'i^B, 
-i{dx + ky)lpA' = £'i>B' - ipA, 



(Al) 



110 

a/e' —aje' —iky/e' —iky/e' 

11 

^iky/e' ~iky/e' p/e' -/3/e' 



(A7) 



The spectrum and the corresponding eigenstates can 



The matrix M is unchanged under the transformation TZ 
and the new "^(x) fulfils the same boundary conditions 
as the old one. 



Appendix B: The transfer matrix 

We denote the wave function to the left of, inside, and 
to the right of the barrier by ^j{x) = GjMjAj, with j = 
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1, 2, and 3, respectively. Further, we have Qi = Qt, and 
Ail =7^3. The continuity of the wave function at a; = 
and X = Wb gives the boundary conditions V'i(O) = V'2(0) 
and 'ip2(Wb) — ''jJsiWb). In explicit matrix notation this 
gives GiAi = ^2^2 and g2M2{Wb)A2 = MiiWb)giA3, 
where Ai = gi'g2M2\Wb)g2'SiMi{Wb)A3. Then 
the transfer matrix J\f can be written as Af — 
gi^g2M2^iWb)g2^giMi{Wb). Let us define A/"' = 
g2M2\Wb)g2\ which leads to Vi(0) = Af'MWb)- 

To treat the case of a (5-function barrier we take the 
limits V ^ oo and Wb such that the dimensionless 
potential strength P = VWb/hvp is kept constant. Then 
CJ2 and A42{Wb) simplify to 



\ 





-iP j 



(Bl) 



/I 


1 








-1 


1 














1 




Vo 





-1 





M2{Wb) 



and M' becomes 






VO 



-iP 






(B2) 



suits; below we give a few examples. In a slightly mod- 
ified notation pertinent to the 2x2 Hamiltonian we set 

a = [— e + ky]^/^, (3 = [e + ky]^/"^, and use the same 
dimensionless units as before. 

Bound states for a single (5-function barrier u{x) = 
PS{x), without accompanying propagating states, are 
possible if fcj, = or fc^ > |e|. In the former case the 
single solution is e = —sign{P)P'^/A. In the latter one 
the dispersion relation is 



e'iP + 2a){P - 2p) + 2P^kl{aP - k^) = 0. (CI) 



The dispersion relation for the KP model obtained 
from the 2x2 Hamiltonian is 



cos(2fcL) -I- 2Fi cos{kL) + F2 = 0, 



(C2) 



where 



' cos P i sin P 

i sin P cos P 

cosP isinP 

i sin P cos P , 



(B3) 



Appendix C: Results for the 2x2 Hamiltonian 

Using the 2x2 Hamiltonian ([3| instead of the 4x4 
one can sometimes lead to unexpectedly different re- 



P P 

Fi = — cosh(/3L) — cosh(Q;L) H sinh(/3L) sinh(Q;L), 

2/3 2a 

F2^^{aPis^ + klPy,) 

+ P cosh(/3L) [a{2e^ - k^) cosh(aL) + e^P sinh(aL)] 

sinh(/3L) [a{e^ - k^/2)P smh{aL) + 2e^ a cosh{aL)] 

(C3) 
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